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Abstract: To determine the magnetic field distribution in the KATRIN main- spectrometer with 
magnetic field sensors that are placed outside the main- spectrometer vessel one can utilize the 
absence of magnetic rotation in main- spectrometer volume. There a scalar magnetic potential V{x) 
can be defined that fulfills the Laplace equation. Large numbers of magnetic field values on an 
outer surface of the main- spectrometer can be sampled by moving and fixed magnetic field sensors. 
These surface samples are used as boundary values in the relaxation of the Laplace equation for 
V{x) and the magnetic field components in the volume. In a simulation involving the KATRIN 
reference solenoid chain, a global magnetic field and an external perturbing solenoid it is shown 
that with this method the original field can be reconstructed within 2 %. 
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systems, architecture, hardware, algorithms, databases). 
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1. The KATRIN setup 

The KArlsruhe TRItium Neutrino experiment flU] (see Fig.[T|) is set up at the Karlsruher Institute 
of Technology (KIT), Germany. It is designed to measure the mass of the electron anti neutrino 
in a direct and model-independent way with a sensitivity of my = 0.2 eV/c^ (90% confidence 
level) from tritium j8 decay flU]- KATRIN uses a magnetic transport field that connects the source 
and detector in combination with integrating electrostatic energy filters (MAC-E-spectrometers). 
Conceptual essentials of the MAC-E spectrometer[0, ED are the magnetic field gradients in pre - 
and main- spectrometer that adiabatically convert cyclotron energy Ecyc into energy Ep parallel to 
the magnetic field lines and vice versa. 



main-spectrometer & EMCS, LFCS 




Figure 1. Schematic view of the KATRIN experiment (total length 70 m) consisting of calibration and 
monitor rear system, with the windowless gaseous T2-source (WGTS), differential pumping (DPS) and cryo- 
trapping section (SPS), the small pre- spectrometer and the large main spectrometer with the large magnetic 
coil systems to compensate the earth magnetic field (EMCS) and to shape the magnetic transport flux (LFCS) 
and lastly the segmented PIN-diode detector. 

At the center of the main- spectrometer (MS) in the minimal magnetic field Ba ^ 3 — 6 jiT, a 
retarding electric field allows an integral energy analysis of Bp. The magnetic field in the analyzing 
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volume defines the magnetic resolution, i.e. the amount of residual cyclotron energy Ecyc that 
can not be analyzed and thus strongly influences the resolution function. Error analysis [M] of 
the influence of uncertainty of the magnetic field in the analyzing plane on the uncertainty of the 
neutrino mass square Am^ leads to a relative accuracy of the magnetic field of ^ < 2.4%. In 
addition, the alignment of magnetic field lines plays a crucial role in the production of secondary 
electrons and electronic background either through penning traps or inner wall contact. 

Large coil systems [ED are arranged around the MS for a) global magnetic field compensation, 
e.g. earth magnetic field (EMCS) and b) fine tuning of the magnetic transport flux with a set 
of large circular low field coils (LFCS) mounted coaxially with the MS (see Fig.|T|). However, 
possible influences of residual external dipoles, magnetization in the MS environment by the high 
field solenoids and/or EMCS, LFCS and the correct orientation of the spectrometer solenoids have 
to be controlled. Due to the extreme MS vacuum conditions the installation of magnetic sensors 
inside the MS is not possible. 



main spectrometer LFCS support ring 
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Figure 2. View of the main spectrometer tank with the LFCS ring system. Right: The mobile sensor unit 
with 2 sensors on the inner belt of a LFCS support ring. 

We therefore propose to determine the magnetic field inside the main spectrometer by taking 
magnetic field samples at an outer surface of the main spectrometer. The sensor network will 
involve fixed position magnetic sensors and mobile magnetic field sensors [0, EL ED which move 
along the inner belts of the LFCS support structure (see Fig.0), close to the outer MS surface but 
well inside the current lines of the EMCS and LFCS . The magnetic field samples serve as boundary 
values for the relaxation of the Laplace equation of the scalar magnetic potential V{x) at the interior 
of the KATRIN main spectrometer. 



2. Volume and surface considerations 

For a volume G with surface area F Amperes equation 

VxB = A,o-(/+^) (2.1) 

can be simplified to the rotationally free case if the current density J is vanishing (7 = 0) and the 
electric field E is constant {dE/dt = 0). 

V X = (2.2) 
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For the KATRIN MS the relevant surface F (see Fig.^ has to be outside the outer MS surface 
and inside the current leading elements (LFCS, EMCS, spectrometer solenoids). As the analyzing 
potential distribution U {x^y^z) inside the MS volume is constant during KATRIN runtime intervals 
(and magnetic field sampling time intervals) the electrical fields produced are time independent. 
Therefore eq. (P^ can assumed to be valid for the KATRIN MS interior. 
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Figure 3. View of the KATRIN main reference solenoid chain and the main spectrometer area. The cylindri- 
cal volume G enclosed by the boundary F. The origin of the coordinate system is located at the geometrical 
center of the MS and the symmetry point of the LFCS. To ensure that all current leading elements are out- 
side the boundary F in y— and z— direction the allowed radius Rq to G is 5600 mm < 7?g < 6155 mm. The 
extreme x— values are -1 1600 mm = —Xn <x< +x„ = 1 1600 mm. 

Vector analysis flTDII states for a scalar function V{x) that: V x V • y (x) = and one can identify 
V{x) with the magnetic scalar potential. 

Utilizing Gauss's law for magnetism V • 5 = we can write down the Laplace-equation (LPE) for 
V{x) 

vV(x,3;,z) = (2.3) 

The finite difference method (FDM) flTTII is chosen to solve the above equation on a 3 dimensional 
rectangular grid, because of its well known numerical stability and the manageable coding effort. In 
the simulation the magnetic field components at a the boundary representing the normal derivatives 
dV /dx^ Bjc; dV/dy = By-, dV/dz = at F can be exported and used in the FD-relaxation as a 
von Neumann boundary values. 



3. Simulation 

The usability of the numerical approach is demonstrated in a simulation based on magnetic field 
values provided by the simulation package PartOpt [0]. The definition of a magnetic scenario (Fig. 
^ at the KATRIN main spectrometer includes: a) the energized KATRIN reference solenoid chain, 
b) the energized LFCS as listed in flT3l1 . c) a magnetic field over G with B^ = 210 mG, By = 35 
mG, B^ = 0, d) a small disturbing magnetic dipole with central induction Be = 600 G adjacent to 
the main spectrometer. 

The field values B^^^By.B^ along the cylindrical surface of volume G with radius Rq ^ 6 m 
between Xmin = —7.03 m < x < Xmax = 6.83 m to cover the cylindrical part of the MS are exported 
in ASCII format. The spacing of the samples in x-direction is 0.45 m in agreement with the real 
x-spacing of the sensor positions. In azimuthal direction a 3° spacing was chosen to get 120 x 2 
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external dipole |B| = 600 G 




Figure 4. PartOpt view of the Simulation scenario. In volume G the effective magnetic field is composed 
of the KATRIN solenoid field, LFCS, an external global field and disturbing external dipole. The perturbed 
magnetic field lines magnetic field lines have been tracked starting from the center of the WGTS. The 
extreme field lines indicate the boundary of the 191 T cm^ nominal magnetic transport flux connecting 
source and detector. 

samples (because 2 sensors are on board) in 15 minutes, the time for one revolution. To simulate 
sensor error the exported values are randomized according to a Gaussean distribution with a 2% 
relative uncertainty. This value was chosen as an upper limit according to the sensor types used in 
[0] . Due to the cylindrical geometry the surface samples points usually do not coincide with surface 
mesh points (cut surfaces problem). Therefore the magnetic samples are interpolated to produce 
values at the regular surface mesh points. The relaxation is performed via a basic 7 point stencil. 
The resulting values for the scalar potential and the values for the magnetic field components are 
generated by deriving V{x) numerically. 
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Figure 5. View of the mesh point structure. Left: G in an interval in , —Xn < x < Xn, Right: on the surface 
1 5 at X — Xj2 •) Xj2 

The relaxation code is written in C. Typically 1400 iterations in 5 minutes on a standard PC 
are performed to meet the terminating condition that the difference for y(0,0,0), the magnetic 
potential at the origin, between successive iterations is < 0.0002. 
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4. Simulation results 



The results of the simulation is displayed as magnetic field components in geometric planes with 
given coordinates within the main spectrometer. Figs.: |5|,|8|,[T(1| show the original PartOpt magnetic 
Borg and the reconstructed magnetic field Brec components for a randomly chosen x^y plane at 
z = 2.4994 m. The relative differences AB are displayed in Fig.: 0,0,|n|,. 




Figure 6. Left: The original magnetic field component Bj^^^ in a in a x,}^ plane at z = 2.4994 m. Right: The 
reconstructed magnetic field values B^^^^ in the same plane. 
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Figure 7. Relative difference between the original B^^^^ and reconstructed B^^^^ magnetic field component 
ABx = {Bx^rg — Bxr^c)/^xorg in a X,}; plane with z = 2.4994 m. The sharp peaks at x ^ 6 arise numerically 
from a division by zero as B^^^^ ^Oin the vicinity of the negatively charged LFCS coil towards the detector 
side as given in [1T3I1 . Elsewhere the difference is less than 2%. 
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Figure 8. Left: The original magnetic field component By^^^ in a in a x,}^ plane at z = 2.4994 m. Right: The 
reconstructed magnetic field values By^^^ in the same plane. The sawtooth structure at the extreme y- values 
are due to the close proximity of the energized LFCS Coils. 




x-Axis [m] 



y-Axis [m] 



Figure 9. Relative difference between the original By^^^ and reconstructed By^^^ magnetic field component 



^y = {By„,-By^J/B. 
division by zero as By^^^ 



y^^g in a X, J plane with z 
0. 



2,4994 m. The sharp peaks arise numerically from a 



-6- 




Figure 10. Left: The original magnetic field component B^^^^ in a in a x, j plane at z = 2.4994 m. Right: The 
reconstructed magnetic field values B^^^^ in the same plane. The sawtooth structure at the extreme y- values 
are due to the close proximity of the energized LFCS Coils. 




Figure 11. Relative difference between the original B^^^^ and reconstructed B^^^^ magnetic field component 
ABz = {Bz^^g ~ ^zrec) l^zorg ^ ^-iJ planc with z = 2,4994 m. The sharp peaks at extreme y -values arise 
numerically from a division by zero as B^ ^ 0. 



Results with similar precision can be found in all areas of the inner volume.Fig. |T^ shows the 
the relative difference AS^ for for a 3;,z-plane at x = 2.475 m. 




Figure 12. Relative difference between the original B^^^^ and reconstructed B^^^^ magnetic field component 
= {^xorg ~^xrec)/^xorg 1^ a plane with X = 2.475 m. The sawtooth structure at the fringes is due to 
the vanishing B — x-component at large radii. 

5. Summary and Outlook 

In a simulation it is shown that with a large number of magnetic field samples taken close to the 
KATRIN main- spectrometer surface and inside the current leading elements of the LFSC -, EMCS 
system and spectrometer solenoids it is possible to determine the magnetic field profile inside the 
spectrometer at least within a 2% precision. With better numerical techniques (e.g. stencils involv- 
ing more meshpoints, interpolation routines with more supporting points) and longer computer 
relaxation times an increase in precission is possible. Also the number and distribution of the sam- 
pling positions on the surface can in the case of the mobile sensor units be varied to achieve better 
results. 

As the front face (at —Xn) and the end face (at of the cylindric volume still intersect the KA- 
TRIN MS volume no samples can be taken there. However, the magnetic field of these surfaces is 
predominantly given by the spectrometer solenoids which can be modeled numerically to produce 
calculated field values. These models can be controlled by fixed position magnetic field sensors 
close to the relevant surfaces. 
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Unlike in a simulation, where the magnetic field components are per se given according to the 
chosen coordinate system, the magnetic field sensors in KATRIN environment have to be aligned 
according to the KATRIN global coordinate system. In the case of moving sensor units moving on 
the inner rails of the LFCS structure as proposed in [0] this requires information about position and 
inclination along the track. 
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